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Abstract 

This article considers the classification of matrix superpotentials that cor- 
responds to exactly solvable systems of Schrodinger equations. Superpotentials 
of the following form are considered: Wk = kQ+P+^R, where k — parameter, 
P, Q and R — hermitian matrices, that depend on a variable x. The list of 
three-dimensional matrix superpotentials is presented explicitly. 
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1 Introduction 



Supersymmetric quantum mechanics presents a powerful tool for finding exact solu- 
tions to problems described by Schrodinger equations [T] . Discrete symmetry between 
a Hamiltonian and its superpartner, called sliape-invariance [2] helps solve the spec- 
tral problem by means of algebraic methods. 

Unfortunately, the class of known problems that satisfy the shape-invariance con- 
dition is rather limited. However, it almost entirely covers the cases when the corre- 
sponding Schrodinger equation can be exactly integrable and has an explicitly pre- 
sented potential. 

A well known paper by Cooper et al [3] presents a classification of scalar potentials 
that correspond to exactly solvable Schrodinger equations. Matrix potentials arise in 
many physics problems. For instance, the motion of neutral nonrelativistic fermion 
that interacts anomaly with the magnetic field generated by a thin current-carrying 
wire, is described by a model with a matrix potential. This model was proposed by 
Pron'ko and Stroganov [1] . Other good examples of problems with matrix potentials 
can be found in works [5]- [7] that describe crystal structures in the Gross-Neve model. 
Some two-dimensional matrix potentials can be found in [8], [9]. Particular cases of 
matrix potentials can be also found in works by [ID]- [IS]- In [11] a particular class of 
two-dimensional matrix potentials is presented. 

A systematic study of the problem was introduced in [15], and [16], where pre- 
sented a classification of superpotentials of the form Wk = kQ + P + ^R, where k — 
parameter, P, Q and R — Hermitian matrices and one of the following conditions is 
held: Q is proportional to a unit matrix, or R is equal to zero. In [T^ a wide class of 
matrix superpotentials was considered, however a full classification was not made. 

This paper considers three-dimensional superpotentials and presents their classi- 
fication. 

2 Solving the spectral problem 

Let us consider the following spectral problem 



where Hk — Hamiltonian with a matrix potential, Ek and ■0 — its eigenvalues and 
eigenfunctions correspondingly. In the Schrodinger equation Hamiltonian has the 
following form 



where Vk{x) — matrix potential dependent on the parameter k and the variable x. 





(2) 
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Let us assume that Hamiltonian can be factorized as follows 



(3) 



where — scalar function of k, that vanishes with a corresponding member in the 
Hamiltonian. From here on, we will drop the sign of unit matrix / and write Ck 
instead of Ckl. As [16] mentions it is enough to consider the operators and a|. of 
the form 



(Ik 



d_ 

dx 



+ Wkix) 



-i^+Wkix), 



(4) 



where Wk — hermitian matrix, that is called a superpotential. 

As Wk is hermitian, operators and a|. are hermitian-conjugate, which allows to 
find the ground state of the spectral problem ([T]), by simply solving the first order 
differential equation. Indeed, multiplying the expression 



alttk-ip = 



(5) 



on the left by the hermitian-conjugate spinor -0^ and integrating the expression on 
the real line M we get 



llflfc^lb = 0, 
where 1 1 ■ 1 1 2 denotes the norm in L2 

ttki^ = 0. 



Hence, 



(6) 



(7) 



The square-integrable function i/^l{x), that is a normahzed solution to the equa- 
tion ([7]), is an eigenfunction of the Hamiltonian, that corresponds to the eigenvalue 

= Cfc, and is called a ground state of the system ([T]). 

Suppose the system satisfies the shape-invariance condition: 



Hi. 



where is a Hamiltonian's superpartner, that is defined by the following formula 



Hi 



ttkol + Ck- 



(9) 



This condition allows to fully discover the spectrum by a series of algebraic operations, 
knowing the eigenvalue of the system i^Kx). Indeed, if we use the condition ([HD, it is 
easy to show that the function 



t„/,o 



x] 



alipt+i{x)\\2 



(10) 
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is a Hamiltonian's eigenfunction with an eigenvalue El = Ck+i- It is called the first 
exited state of the system ([T]). Analogously, by induction we prove that the function 




■ ■ ■ a 



■ ■ ■ a 



t 

'k+n-l 



'k+n-1 




(11) 



is a Hamiltonian's eigenfunction with an eigenvalue = Ck+n- It is called the n 
exited state of the system ([T]). 

Consequently, if the system of Schrodinger equations ([T]) satisfies the shape- 
invariance condition, it can be exactly integrated. 

3 The classification problem 

Since the shape- invariant potentials correspond to the exactly integrable systems 
of Schrodinger equations, it would be interesting to widen the class. Let us find 
all the Hamiltonians that allow factorization (|3]) and satisfy the shape-invariance 
condition (|8]). In terms of the superpotential these conditions can be written through 
a single equation 



where = c^+i — c^. Thus, to solve the given problem, it is enough to find all the 
superpotentials, that satisfy the equation (IT^ . 

In the general form, the problem is rather complex. However, it can be solved if 
we consider superpotentials of a chosen class. This paper considers superpotentials 
of the form 



where P, Q and R — Hermitian matrices of 3 x 3 dimension. Additionally, let us 
consider that Q is not proportional to a unit matrix and R is not a zero matrix. 
These cases were considered in [15], [16]. 

We are interested in the irreducible superpotentials, i.e. the ones that can not 
be reduced to a block-diagonal form by means of an unitary transformation, that 
does not depend on variable x. In the reducible superpotentials case, the problem is 
divided into a set of similar problems of a smaller dimension. 

In the following section the equation describing unknown matrices P, Q and R 
is obtained and solved, given the corresponding superpotentials satisfy the equa- 
tion ([T^. 



wl + wl = wl^,-wu, + Q 



(12) 



Wk = kQ + P + -R, 



(13) 
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4 The determining equations 

To get a system of determining equations, let us substitute the expression f[T^ . for 
the superpotential, into the equation (IT^ and separate the variables. Hence, 



R' = 0, 



{P, R} + x = 0, 



X {2k + l)u^ 



+ 



k{k + l) A;2(A; + 1)2 



(14) 
(15) 

(16) 
(17) 
(18) 

(19) 



where u, cu, x — arbitrary real constants. 

As previously shown in [16|, matrix Q can be diagonalized by means of a unitary 
transformation that does not depend on x. Then the expression (HM is reduced to a 
system 



i + 



1...3, 



(20) 



where qi — diagonal elements of the matrix Q. 
It has the following solutions: 

gj = Atan(Ax + 7j), z = 1...3, z/ = A^; 
— Atanh(Aa; + 7j), i = l...m 
— A coth(Ax + 7j), i = m + l...l 
±A, i = l + !...?> 

^ " — 1 . . .m 

m + 1 . . . 3 



(21) 



0, i 



z/ = 0, 



where 7^ G M, i = 1 . . . 3 — integral constants, m, / = . . . 4, gang a . . .b, where a > b 
is considered empty. 

As matrix Q is diagonal, then the linear equation (|T5|) is dissolved and can be 
solved element-wise: 



Pa = J tan(Ax + 7^) + (pu sec(Ax + 7^), z = 1 . . . 3 

Pij = fij \J sec(Ax + 7j) sec(Ax + 7^-), i = 1 . . .?>, j = 1 . . .?> 



{22) 
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• liu = -A^ 



Pa 



Pij 



— ^ tanh(Ax + 7i) + sech(Ax + 7j), i 
coth(Ax + 7i) + v^ii csch(Ax + 7i), i 

i = 1 
i = 1 
i = 1 



: 1 . . . m 

■■ m + 1 . . .1 

:/ + !.. .3 

. m, J = 1 . . . m 
.m, j = m + 1 . . .1 
.m, j = I + 1 . . .3 



(fij ^ysech{Xx + 7i) sech(A2; + 7j), 
(/jjj A/sech(Ax + 7j) csch(Ax + 7^), 
ipij ^^/sech{Xx + 7j) exp(±Ax), 
(/jjj A/csch(Ax + 7j) csch(Ax + 7^), i = m + l...l,j = m + l...l 
(fiij \/ csch(Ax + 7j) exp(±Ax), 
(/)yexp(±Ax), {-k) 

^ij, (*) 



m+l.../,j = / + !.. .3 
/ + 1...3,i = / + 1...3 
/ + 1...3,i = / + 1...3 



(23) 



Let us use the (*) formula, when the corresponding diagonal elements qi and qj have 
same signs, and (*), if the signs are different. 

• If z/ = 



Pi', 



Pij 



IJ.X 



X + 27,) 



x + % 

-fiX + ifii, 

Vij 

^(a; + 7i)(x + 7j)' 
'fij 



y/xTYi'' 
Vij, 



i = 1 . . .m 
i = m + 1 . . .3 
i = 1 . . .m, j = 1 . . .m 

i = 1 . . .m, j = m + 1 . . .3 
i = m + 1...3,j = m + 1...3 



(24) 



where tpji = Tp~j e C — integral constants, and m and / in the gangs correspond to 
such in the formula (l^T]) . 

Through equations (ITBil . and (ITTI) we can conclude that i? = (rj^) is a constant 
matrix, square of which is proportional to a unit matrix. 

The last equation f|T8|) imposes extra conditions onto the constants /i, x and in- 
tegral constants (fij,rij. Let us show that, in the cases when Q is not a constant 
matrix, the following condition has to be held 



H = 0, X = 0. 



(25) 
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Let us consider the elements {P, R}ij in the equation f lTS]) . that correspond to 
a nonconstant element of the matrix Q: 



4:^riSi{x) + '^{vipipip + rip(pip)r]ip{x) = -x, 

3 

'^■imMii.x) + ^j(x)) + ^(^rjp?7ip(x) + ipjpr-r]jp{x)) =0, j 



(26) 



1...3, 



where rjij 



p=l 

multipliers of cpij in the matrix P; 

A2 



1 ... 3 defined by the formula 



i^tan(Ax + 7i), 

— J tanh(Ax + 7i 

— J coth(Ax + 7j~ 

2{x + 70 ' 



V 
V 



-A2 
-A2 





(27) 



and similar to ^j, or the constant. 

Since ^i{x),'r]ij{x) and 1 are linearly independent with each other, then the sys- 
tem 0261) will be consistent if only fi and x are equal to zero, or the entire column rjj 
of the matrix R is equal to zero. However, in the last case matrix i? is a singular 
matrix, which contradicts the condition f|T7|) . hence fl25|l is proven. 

The achieved conditions significantly simplify the equation f lTSj) and allow to solve 
it element-wise, simplifying the results by means of unitary transformations that do 
not depend on x. 

In the following section the solutions to the ffTSl) equation and results (1^ — fl2^ 
are collected and presented in the form of a list of the three-dimensional superpoten- 
tials. 



5 Three-dimensional matrix superpotentials 

Let us write down matrices P, Q and R separately for convenience. For matrices P 
and R basis of the Gell-Man type is used: 



ei 



64 



^0 1 0' 

1 I ee 

^0 0, 

^0 r 

) eg 

1 0, 




^0 -I 0' 

I I 65 

.0 0^ 

'1 0^ 

-1 ) 67 

,0 0. 




(28) 
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and Co defines unit matrix /. And for matrix Q we use a diagonal form Q ~ 
diag {gi, (72, gs} = + ^1^63 + ^^tt^eg. 

The list of nonequivalent representations of the matrix Q: 

• liu = -A2 

Q = diag {Atan(Aa; + 71), Atan(Aa; + 72), Atan(Ax + 73)} (29) 

• If = A^ 

Q — diag{— Atanh(Aa; + 71), — Atanh(Ax + 72), — Atanh(Ax + 73)} (30) 

Q — diag{— Acoth(Aa; + 71), — Acoth(Aa: + 72), — Acoth(Ax + 73)} (31) 

Q — diag{— Atanh(Aa; + 71), — Atanh(Ax + 72), — Acoth(Ax + 73)} (32) 

Q = diag {—A coth(Ax + 71), —A coth(A.T + 72), — Atanh(Aa; + 73)} (33) 

Q — diag{— Atanh(Aa; + 71), — Atanh(Ax + 72), ±A} (34) 

Q — diag{— Acoth(Aa; + 71), — Acoth(Aa: + 72), ±A} (35) 

Q — diag {± A, — Atanh(Ax + 72), — Acoth(Ax + 73)} (36) 

Q = diag {±A, ±A, — Atanh(Ax + 73)} (37) 

Q = diag{±A,±A, -Acoth(Aa; + 73)} (38) 

g = diag{A, -A,-Atanh(Ax + 73)} (39) 

Q = diag{A, -A, -Acoth(Ax + 73)} (40) 

• If = 

Q = diag I I (41) 

1^ a; + 71 X + 72 a; + 73 J 

g = diag( ^,o| (42) 

[ a; + 7i a; + 72 J 

Q = diag|o,0, ^1 (43) 

I a; + 73 J 

Where two cases are considered: 71,72,73 — all are different and 71 = 72 7^ 73- 
The list of nonequivalent representations of the matrix R for different 71,72,73: 

R = ±oj(e3 + eg) (44) 

it! = £62 + res i (45) 
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R = ±u{e3 - eg) (46) 

R=^{r± u)eo + ^irT uj)e3 + ee^ - ^(3r ± co)e8 (47) 

R = ±uj{eo - 268) (48) 

R=^{r± u)eo - ^{r t w)e3 + £67 - ^(3r ± u)e8 (49) 

and for 71 = 72 7^ 73: 

R = ±u{ei ± eg) (50) 

11 1 
R = -{r± uj)eo + ^{r T i^)e3 + ee^ - -(3r ± u)e8 (51) 

where p,r,(p,e — real constants, p ^ and + = w^. In the f l5Up case, signs 
before ci and eg are chosen independently. 

Below the list of nonequivalent representations of the matrix P is presented. 

For different 71,72,73, matrix R, that is defined by formula i^^, corresponding 
to Q, that is defined by formulas ([2n])-(I3S]), (EHD-dll) 

P = 4)61^/ sec(A2; + 71) sec(A2; + 72) + pe^^/ sec(Ax + 72) sec(Ax + 73) (52) 



P = 4>e\\/ sech(Ax + 71) sech(Ax + 72) + pe^ sech(Ax + 72) sech(Ax + 73) (53) 
P = 4>e\\/ csch(Aa; + 71) csch(Aa; + 72) + pes\/ csch(Aa; + 72) csch(Ax + 73) (54) 



P = 0eiv^ sech(Ax + 71) sech(Ax + 72) + pee\/ sech(Ax + 72) csch(Ax + 73) (55) 
P = (pei^y csch(Aa; + 71) csch(Ax + 72) + pcq-s/ csch(Ax + 72) sech(Aa; + 73) (56) 



P = 4>ei^/ sech(Ax + 71) sech(Ax + 72) + peQ^y sech(Ax + 72) exp(±Aa;) (57) 
P = (pei-s/ csch(A2; + 71) csch(Aa; + 72) + peQ\/ csch(Aa; + 72) exp(±Ax) (58) 



P = 4>ei^/ exp(±Ax) sech(Ax + 72) + pcg^/ sech(Ax + 72) csch(Ax + 73) (59) 

P = (pCi + pcQs/ sech(— Ax) sech(Aa; + 73) (60) 

P = 0ei + peQ\/ sech(— Ax) csch(Aa; + 73) (61) 

P= ^ "^'^ + , (62) 

V(x + 7i)(x + 72) a/(x + 72)(x + 73) 

p= , "^'^ +-m= (63) 

v/(x + 7i)(x + 72) V^+l2 
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For different 71,72,73, matrix R, that are defined by formula (jH]), corresponding 
to Q, tfiat is defined by formulas ([2n])-(I3S]), (EHD-dll) 



P = pei\/ sec(Ax + 71) sec(Ax + 72) 


(64) 


P = pe\\/ sech(Ax + 71) sech(Ax + 72) 


(65) 


P = V^\\/ csch(Ax + 7i) csch(Ax + 72) 


(66) 


P = "pew/ secli(Ax + 71) sech(Ax + 72) 


(67) 


P = peis/ csch(Aa; + 71) csch(Ax + 72) 


(68) 


P = pe\\/ secli(Ax + 71) sech(Ax + 72) 


(69) 


P = pe-\\/ csch(Ax + 7i) csch(Ax + 72) 


(70) 


P = TiP^ \ / P'^wi -\~\t \ ^pplif At -I- ^ol 




P = pci 


(72) 


P = pei 


(73) 


v/(x + 7i)(a; + 72) 


(74) 


p_ Pei 

v/(x + 7i)(x + 72) 


(75) 



For different 71,72,73, matrix R, that is defined by formula fHBl) . corresponding 



to Q, that is defined by formulas (I5UD. dlU]) 

P = 4>ei^/ exp(±Ax) sech(Ax + 72) + pe4-\/ exp(±Ax) csch(Aa; + 73) (76) 

P = 0ei + pe4^/ exp(Ax) sech(Ax + 73) (77) 

P = (pci + pe^^J exp(Aa:) csch(Ax + 73) (78) 

For different 71,72,73, matrix i?, that is defined by formula (jUj), corresponding 
to g, that is defined by formulas ([S2D-(I3nD, (15^]). (|iU|), 

P = ptt^ sech(Aa; + 71) csch(Ax + 73) (79) 

P = pe/^\J csch(Ax + 7i) sech(Ax + 73) (80) 

P = pe^s/ sech(Ax + 71) exp(±Ax) (81) 

P = pe/^\J csch(Ax + 7i) exp(Ax) (82) 
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P = 'pe/^\J exp(±Ax) csch(Aa; + 73) (83) 



P = "pt/^^J exp(Ax) sech(Ax + 73) (84) 
P = pe^-s/ exp(Ax) csch(Ax + 73) (85) 
P = (86) 

For different 71,72,73, matrix R, that is defined by formula fHSl) . corresponding 
to Q, that is defined by formulas ([S2D-(I3nD, (15^]). dlCT]). (112]) 

P = 4>e4\/ sech(A2; + 71) csch(Ax + 73) + pe^^/ sech(Aa; + 72) csch(Ax + 73) (87) 



P = 4>e4\/ csch(Ax + 7i) sech(Ax + 73) + peQ^y csch(Ax + 72) sech(Ax + 73) (88) 
P = (^e^ sech(A2; + 71) exp(±As) + pe^^J sech(Ax + 72) exp(±Aa;) (89) 
P = (^e^ csch(A2; + 71) exp(A2:) + pe^\J csch(Aa; + 72) exp(±Aa;) (90) 



P = 4>e4^/ exp(±Ax) csch(Ax + 73) + pe^^/ sech(Ax + 72) csch(Ax + 73) (91) 
P = 4>e4\/ exp(Ax) sech(Ax + 73) + peQ\/ exp(— Ax) sech(Ax + 73) (92) 



P = 4>e4^/ exp(A2;) csch(Ax + 73) + pe^^/ exp(— Ax) csch(Ax + 73) (93) 

P= ^ +P^ (94) 

For different 71,72,73, matrix R, that is defined by formula fH^ . corresponding 
to Q, that is defined by formulas (15^]). dlU]) 

P = peQ-s/ exp(±Ax) csch(Ax + 73) (95) 



P = V^&J exp(Ax) sech(Ax + 73) (96) 

P = peQ-s/ exp(Ax) csch(Ax + 73) (97) 

For 7i = 72 7^ 73, matrix R, that is defined by formula (150|) . corresponding to Q, 
that is defined by formulas d^-dSSD, §7^, ([381), (ill)-(!l3l), where the sign in the 
matrix P before eg — opposite to the chosen sign in the matrix R before eg 



P = pes sec(Ax + 71) + ^(64 =F ee)^ sec(Ax + 71) sec(Ax + 73) (98) 



P = pes sech(Ax + 71) + ^(64 =F 5^)^/ sech(Ax + 71) sech(Ax + 73) (99) 
P = pes csch(Ax + 71) + ^(64 =F ee)^/ csch(Ax + 71) csch(Ax + 73) (100) 



P = pes sech(Ax + 71) + ^(64 =F ee)^/ sech(Ax + 71) csch(Ax + 73) (101) 
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p 


= ne-x cschfAx + 71) + MeA =F eeK/ cschfAa; + 71) sechfAa; + 73) 


(102) 


p 


= pes sechfAx + 71) + 0(64 =F eelx/ sechfAx + 71) exp(±Ax) 


(103) 

V / 


p 


= pes csch(Ax + 71) + ^(64 =F e6)A/ csch(Ax + 71) exp(±Ax) 


(104) 


p 


= pea exp(±Ax) + ^(64 =F e6)A/ exp(±Ax) sech(Ax + 73) 


(105) 


p 


= pes exp(±Ax) + 0(e4 =F ^&)\/ exp(±Ax) csch(Ax + 73) 


(106) 


p 


_ pes 1 0(e4=Fe6) 
2^ + 71 V + + 73) 


(107) 


p 


_ peg ^ 0(64 =F eg) 
2^ + 71 + 7i 


(108) 


p 


0(e4 =F ee) 
= pes + ^ 


(109) 



For 7i = 72 7^ 73, matrix i?, that is defined by formula flSTj) . corresponding to Q, 
that is defined by formulas ([2n])-(I3SD, dSZD. dM]), (gT])-(!i51) 



p 


= pe^^J sec(Ax + 71) sec(Ax + 73) 


(110) 


p 


= pei\/ sech(Ax + 71) sech(Ax + 7s) 


(111) 


p 


= pe/^^J csch(Ax + 71) csch(Ax + 7s) 


(112) 


p 


= pei\/ sech(Ax + 71) csch(Ax + 7s) 


(113) 


p 


= pe^^J csch(Ax + 7i) sech(Ax + 7s) 


(114) 


p 


= pe4v^ sech(Ax + 71) exp(±Ax) 


(115) 


p 


= pe^is/ csch(Ax + 71) exp(±Ax) 


(116) 


p 


= pei\/ exp(±Ax) sech(Ax + 7s) 


(117) 


p 


= pe^^J exp(±Ax) csch(Ax + 7s) 


(118) 


p 


pe4 

v/(x + 7i)(x + 73) 


(119) 


p 


_ pe4 
a/^ + 7i 


(120) 


p 


pe4 
y/x + 7s 


(121) 
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The corresponding three-dimensional matrix superpotentials can be easily recre- 
ated by means of the formula f|T3|l . 

6 Conclusions 

The article intended to find three-dimensional matrix superpotentials, that corre- 
spond to exactly integrable systems of Schrodinger equations. In the general form 
the problem remains rather complex to analyze, but it was solved by limiting the 
class of the considered superpotentials by a simple case f|T3|) . 

The list of considered three-dimensional superpotentials is presented by means of 
matrices Q (formulas ([22]) -(USD), R (formulas dSD-dSI])) and P (formulas 
Corresponding superpotentials is easy to get explicitly by means of formula (|T3l) 
beginning with formulas for the matrix P. 

The received results can be interesting for physicists, in order to find new exactly 
integrable physical models. 

Even though the problem was solved entirely, the achieved results are limited 
by the form of superpotential f|T3|) and the dimension of the superpotential. In the 
upcoming works, a more general class of superpotentials will be considered. 
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